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1.  Introduction 


One  of  the  earliest  applications  of  queueing  theory  to  spares 
provisioning  problems  was  the  work  of  Taylor  and  Jackson  [121,  in  which 
a finite  source  queueing  model  with  spares  was  used  to  determine  the 
number  of  spare  engines  required  to  maintain  a fleet  of  aircraft  at  a 
certain  efficiency  level.  A bibliography  of  some  recent  work  on  queue- 
ing approaches  to  provisioning  type  problems  is  given  in  bureau  [6].  Of 
particular  interest  are  [2],  [8],  and  [9].  The  problem  treated  in  this 
paper  is  the  determination  of  an  adequate  number  of  spares  and  repair  lines 
(servers)  for  replacing  and  repairing  components  which  randomly  fail,  as- 
suming that  the  failed  components  are  replaced  by  spares  (if  available)  and 
once  repaired,  in  turn  become  spares.  A multi-year  planning  horizon  is  con- 
sidered, allowing  for  growth  both  in  component  population  size  and  component 
reliability. 

When  a component  fails,  a request  for  a spare  is  immediately  placed. 
If  no  spare  is  available,  a delay  occurs.  A service  level  constraint  of 
10%  is  imposed  on  such  delays;  that  is,  we  desire  a capability  such  that  at 
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at  least  90%  of  the  requests  for  spares  are  immediately  filled  from  on- 
hand  spares  inventory  (at  most  10%  backordering  of  spares  requests).  This 
service  criterion  is  often  referred  to  as  availability  or  fill  rate. 

The  objective  then  becomes  one  of  minimizing  expenditures  for  spares 
and  servers  subject  to  a 90%  fill  rate  constraint.  The  classical  machine 
repair  queueing  model  with  spares  is  the  driving  model  in  that  it  determines 
fill  rate  given  a specific  number  of  spares  and  servers,  the  component  popu- 
lation size,  failure  rate  and  service  rate.  Times  between  failures  are 
assumed  to  be  exponentially  distributed  random  variables  as  are  service 
times  (which  include  the  time  required  for  removal,  transportation,  and 
repair).  This  queueing  model  is  then  embedded  in  a heuristic  cost  optimi- 
zation model  which  determines  a "good"  mix  of  spares  and  servers  for  each 
year  in  the  planning  horizon,  while  satisfying  the  fill  rate  constraints. 

An  alternate  queueing  model  is  also  presented  which  treats  removal, 
transportation,  and  repair  as  separate  serving  stations,  but  which  must 
assume  an  "infinite"  calling  population.  The  accuracy  of  this  assumption 
is  also  investigated. 

The  methodology  presented  here  has  been  used  for  provisioning  servers 
and  spares  for  a fleet  of  marine  gas  turbine  engine  ships  (each  engine 
having  two  components — a gas  generator  and  a power  turbine).  It  is,  of 
course,  applicable  to  other  similar  provisioning  problems. 

2 . Queue ing  Model 

Tn  order  to  determine  the  probability  of  a request  for  a spare  being 
met  without  delay,  it  is  necessary  to  calculate  the  probability  of  various 
numbers  of  components  in  the  repair  system  at  any  particular  time.  Changes 
in  the  population  with  respect  to  size  and  reliability  take  place  at  random 
times  throughout  each  year  so  that  unless  (or  until)  population  size  and 
component  reliability  stop  changing,  steady  state  cannot  be  approached.  It 
appears  analytically  intractable  to  calculate  transient  state  probabilities. 
As  an  approximation,  we  consider  the  population  to  be  in  steady  state  at  its 
average  size  and  reliability  for  an  entire  year,  changing  instantaneously  to 
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a new  steady-state  position  at  new  average  values  at  the  beginning  of  each 
new  year.  In  situations  where  failures  are  frequent,  transient  effects 
should  die  out  quickly  and  our  steady-state  approximations  should  be  ade- 
quate. However,  if  failures  are  infrequent,  the  transient  effects  will 
take  longer  to  disappear  and  might  present  a problem.  Nevertheless,  for 
most  provisioning  purposes,  the  assumption  of  instantaneous  steady  state 
should  be  adequate. 

At  any  point  in  time,  the  population  is  composed  of  units  which 
may  have  different  failure  rates,  since  units  added  to  the  population  in 
later  years  generally  are  more  reliable  due  to  technological  learning. 

It  is  assumed  here  that  for  each  year,  all  components  have  identical  failure 
rates  equal  to  the  population  average,  which  changes  on  a yearly  basis,  as 
described  above.  This  assumption  that  al 1 components  operate  at  the  popu- 
lation average  failure  rate  was  investigated  in  Reference  [5],  and  the 
results  will  be  briefly  described  la^er  in  this  paper. 

We  introduce  the  following  notation: 

p^  ^ = steady-state  Pr{n  components  in  repair,  year  i} 
c^  = number  of  repair  facilities,  year  i 
y^  = number  of  spares,  year  i 

A^  = component  failure  rate  (Poisson  mean),  year  i 

A^  = population  average  failure  rate,  year  i 

R.  = expected  number  of  components  repaired,  year  1 

= component  population  size,  year  i 

1/y.  = average  service  (removal,  transportation,  and 
repair)  time,  (exponential  mean),  year  i. 

The  equations  for  the  steady-state  probability  of  n items  "down" 
are  given  as 
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Equation  Set  (1)  gives  the  standard  formulas  for  a machine  repair  model 
with  spares  (see  Gross  and  Harris  [4],  page  123,  for  example). 

Using  the  average  failure  rate,  A , for  each  year  allows  for  the 

incorporation  of  changing  component  reliability  as  the  years  progress. 
Equation  Set  (2)  shows  how  A is  computed.  It  is,  perhaps,  easiest  to 
explain  in  the  context  of  the  marine  gas  turbine  engine  example  mentioned 
previously.  Year  one  begins  with  the  first  introduction  of  a few  gas  tur- 
bine ships  and  each  succeeding  year  brings  into  the  fleet  additional  ships 
until  the  fleet  reaches  full  strength.  The  components  introduced  in  the 
later  years  generally  have  improved  reliability,  either  through  learning 
or  a conscious  component  improvement  program  (CIP).  Thus  A^  tends  to 

be  smaller  than  A^_^  • When  the  population  size  is  increasing  (N^  > 

j ) , the  average  failure  rate  over  all  components  in  the  population  for 
year  i,  A^  , is  given  as  the  weighted  average  of  the  new  components  in- 
troduced into  the  population  at  the.  best  current  achievable  failure  rate, 

A^  ; those  repaired  during  the  past  year  at  that  year's  best  achievable 

failure  rate,  A^_^  5 and  those  old  components  in  the  population  which  were 

not  repaired  and  are  thus  operating  at  the  old  average  failure  rate, 

A^_^  . If  the  fleet  size  is  decreasing  (N^  < , that  is,  some  ships 

may  be  retired  from  service,  then  the  computation  is  given  by  the  second 
equation  of  Equation  Set  (2),  namely,  the  weighted  average  of  those  repaired 
last  year  coming  in  with  a failure  rate  of  A^  and  those  not  repaired  but 

operating  at  the  old  average  A^_^  • The  average  number  repaired  in  a year 

is  given  by  Equation  (3). 


The  assumption  that  all  components  operate  at  an  average  failure 
rate  was  investigated  in  Reference  [5]  by  developing  the  exact  model  for 
unequal  failure  rates  with  N=1  , Y=1  , c=l  , so  that  one  component  has  a 


Aj  while  the  other  has  a failure  rate  of 


failure  rate  of 
out  that  (i)  assuming  both  units  operate  at 


It  turns 


Xl+X2 


gives  conservative 


5 


T-344 


results  with  respect  to  availability,  i.e.,  the  actual  availability  is 
greater  than  that  computed  using  A and  (ii)  even  for  sizable  differences 
in  A^  and  ■>  t*'e  approximate  availability  using  A is  close  to  tiie 

actual.  The  reader  is  referred  to  Reference  [5]  for  more  detail.  Result 
(i)  has  some  intuitive  appeal  since  components  having  higher  failure  rates 
would  be  expected  to  fail  (and  hence  be  in  repair)  more  often.  Thus,  the 
more  reliable  components  operate  a larger  proportion  of  the  time  so  that 
the  actual  failure  rate  for  the  system  would  tend  to  be  lower  than  the  arith- 
metic average  of  the  failure  rates  over  all  components.  Nevertheless,  it 
appears  from  Result  (ii)  that  this  effect  is  not  highly  significant  and  that 
using  the  population  average  failure  rate  for  each  component  is  an  adequate 
approximation  in  most  cases.  For  example,  using  the  exact  model  with  N = 

Y = c = 1 , it  was  found  that  even  for  A ^ twice  that  of  A^  , the  percent- 
age error  in  availability  was  no  more  than  5%  when  using  A as  the  failure 
rate  for  each  component.  In  actual  applications,  the  individual  component 
failure  rates  tend  to  be  much  closer  together  than  a factor  of  two,  since 
reliability  growth  is  gradual,  and  hence  we  feel  confident  that  the  average 
failure  rate  assumption  is  justifiable. 

Once  the  p . are  obtained,  the  fill  rate  constraint  is  given  by 
n,  l 

(temporarily  dropping  the  subscript  i ) 

v-1 

L > °-90  , 

n=0  n 


since  when  n components  are  down,  there  are  y-n  spares  available 

(n  < y)  . Thus  the  probability  of  no  spares  on  hand,  P , is  1 - 

out 


y'1 

? . p and  the  percentage  of  requests  filled  immediately  from  on-hand 


n=0 

spares  is  A(l-P  )/A  = 1 - P 

out  out 


y-1 

- £ pr 

n=0 


However,  the  probability  that  a failed  component  finds  n in  the 
system  should  be  used  in  the  fill  rate  computation.  In  common  queueing 
terminology  these  are  the  "arriving  customer"  probabilities  and,  in  the 
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context  of  this  paper,  correspond  to  the 
which  generate  requests  for  spares  and 
failure  point  probabilities.  For  the  T 
sidered  in  this  paper  the  failure  point 
general  time  probabilities  given  by  Equ 
no  spares  queue  (see,  e.g. , Cooper  [1], 
abilities  are  equivalent  to  the  genera] 
calling  population  of  one  less  but  this 
spares  case. 

The  failure  point  probabilities 
queue  (denoted  by  q^  as  contrasted  tc 

denoted  by  p ) may  be  derived  as  foil 


; occurrence  of  component  failures 
:hus  they  shall  be  referred  to  as 
Lnite  source  with  spares  queue  con- 
probabilities  are  not  equal  to  the 
it  ion  Set  (1).  For  the  finite  source- 
pp.  82  ff.)  the  failure  point  prob- 
timo  probabilities  for  a finite 
relationship  does  not  hold  for  the 

for  the  finite  source  with  spares 
the  general  time  probabilities, 

ows.  Using  Bayes'  theorem. 
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For  the  R calculation,  the  general  time  probability  is  required 
so  that  Equation  (1)  is  used  as  given,  and  in  fact  R equals  the  denomina- 
tor term  of  q^  multiplied  by  A . 


3.  Cost  Model 

The  cost  model  considers  four  types  of  costs:  (1)  purchase  cost 

of  spares,  (2)  purchase  cost  of  service  channels,  (3)  repair  costs  and 
(4)  investment  costs  in  component  improvement  (reliability  improvement) 
programs.  Annual  operating  costs  associated  with  running  the  spares  in- 
ventory system  or  operating  the  service  channels  are  not  included.  While 
the  variable  portion  of  these  costs  can  be  important,  they  are  not  consid- 
ered in  this  model  since  the  major  purpose  of  such  a strategic  planning 
model  as  this  is  for  capital  budgeting  and  hence  it  is  the  purchase  ex- 
penditures which  are  of  prime  concern  at  this  stage  of  the  decision-making 
process. 

The  optimization  problem  is  an  integer-nonlinear-programming  problem 
with  a hard- to-manage  objective  function  and  highly  complex  constraints. 

An  integer  programming  algorithm  approach  has  not  as  yet  been  successful  in 
solving  the  problem  and  hence  a heuristic  method  is  utilized. 
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The  heuristic  cost  "optimization"  algorithm  considers  explicitly 
only  the  first  two  categories  of  costs,  although  the  repair  and  component 
improvement  program  (CIP)  costs  are  always  calculated  for  each  year  so 
that  sensitivity  analyses  with  respect  to  various  CIP  programs,  service 
times,  etc.,  can  be  performed.  it  should  be  noted  that  even  with  no  CIP 
program,  A and  N would  still  change,  hut  the  CIP  program  Influences 
the  magnitudes  of  these  changes.  The  present  worth  of  the  cost  stream 
over  the  planning  horizon  is  always  calculated  for  just  such  a use.  This 
will  be  illustrated  in  the  next  section. 


The  cost  minimization  problem  can  be  stated  as  follows.  Considering 
the  additional  notation, 

C.  , = purchase  cost  of  a repair  channel,  year  i 

M * 

C„  , = purchase  cost  of  a spare  component,  year  i 

^ 9 

C„  = cost  of  repairing  a component,  year  i 

C.  . = investment  in  reliability  growth  (CIP)  program, 
year  i 

d = discount  factor 

k = number  of  years  in  planning  horizon. 


we  desire  to 


Minimize  Z 

Cl*  ^"2 ’ " " ’ ’ Ck’  yi,y2’ * ' * ,yk 


y-l 


Subject  to 


d [Cia(ci  ~ 

Ci-1)  + C2,i(yi  ~ yi-l) 

■C3,iVC4,i] 

(4) 

q > 0.90  , 

1 n,i 

(i=l,2,...,k)  (5) 

>_  0 ; integer 

0 ; integer. 

The  plus  superscript  on  the  first  two  cost  terms  indicates  the  term  is 
zero  if  the  factor  in  parentheses  is  negative,  that  is. 
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(a-b)+  = max{(a-b),0[  . 

The  last  term  In  Equation  (4),  C,  . , Is  the  cost  of  the  CIP 

4,  i 

and  does  not  explicitly  depend  on  c.  or  y.  , although  indirectly  one 

can  argue  that  if  C.  . is  increased,  X.  will  be  smaller  thus  affecting 

4,1  l 

the  p . calculation  and  hence  the  constraint  [Equation  (5)].  The  heu- 
n , l 

ristic  algorithm  ignores  this.  This  CIP  effect  is  observed  via  a sensi- 
tivity type  analysis  referred  to  above. 


The  third  term,  which  involves  R.  , is  a function  of  y^  [ex- 
plicitly; see  Equation  (3)|,  and  y.  and  c implicitly  through  p . . 

l i n , l 

However,  when  compared  to  purchase  costs  of  spares  and  servers,  the  annual 
repair  costs  should  be  small  and  are  ignored  by  the  heuristic  algorithm. 

It  is  included  in  the  cost  calculation  since  it  is  directly  influenced  by 
reliability  growth  and  is  needed  for  any  sensitivity  study  of  CIP.  Thus 
the  heuristic  algorithm  looks  only  at  the  costs  of  purchasing  spares  and 
servers  and  furthermore,  considers  only  one  year  at  a time,  that  is,  it 
attempts  to 


Min 

Vyi 


zi  = cl,i(ci- 


Ci-1)+  + C2,i(yi  " 


yi-l} 


for  each  i=l,2,...,k  , in  a sequential  manner,  using  the  "best"  combina- 
tion it  finds  for  year  i-1  (c^  ^ ,y^  as  the  starting  point  for  its 
search  for  • 

For  the  year  i computations,  the  algorithm  first  checks  to  see  if 
c i _ 1 , y 1 | satisfy  the  constraint.  If  the  constraint  is  not  met,  the  al- 
gorithm computes  the  ratio  A = C„  {/C,  • and  takes  the  largest  integer 

value  contained  therein.  (In  the  applications  we  have  studied,  C . was 

^ 1 

always  larger  than  C . If  this  is  not  the  case,  the  algorithm  should 

1 » * 

be  modified  to  compute  C /C.  . > and  the  words  "server"  and  "spare" 

L , 1 L , 1 

interchanged  in  the  description  that  follows.)  Thus  for  an  (approximately) 
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r 

equal  dollar  expenditure,  we  can  purchase  one  spare  or  A servers. 

Both  points  (c^_^  + y1_1)  and  (c^_^,  y^  j + 1)  are  decked  to 

determine  the  fill  rate  percentage  via  use  of  the  queueing  model,  and 
the  point  with  the  largest  fill  rate  percentage  is  chosen.  This  pro- 
cedure is  continued  from  the  "new  points"  chosen  until  the  fill  rate 
reaches  (or  exceeds)  90%.  If  90%  is  exceeded  (and  it  generally  will  be 
because  of  the  integer  values  required  for  c and  y)  a reduction  pro- 
cedure is  employed  to  see  if  any  servers  can  be  decreased  (since  A 
servers  at  a time  have  been  added)  and  still  maintain  a 90%  fill  rate. 

Also,  if  the  feasible  region  has  been  entered  by  adding  A servers  (as 
opposed  to  adding  a spare)  in  addition  to  decreasing  servers,  an  attempt 
is  made  to  reduce  the  number  of  spares.  This  reduction  portion  of  the 
algorithm  is  intended  to  find  a solution  as  near  to  the  constraint  bound- 
ary as  possible.  The  procedure  is  schematically  shown  in  Figure  1,  which 
displays  all  the  points  (c.,y)  evaluated  by  the  algorithm  along  with  the 
path  of  greatest  improvement  in  fill  rate.  Note  for  each  point  on  the 
path,  fill  rate  is  calculated  for  the  points  immediately  above  and  to  the 
right,  and  then  a move  is  made  to  the  one  with  the  largest  fill  rate. 

I 

Should  the  starting  value  for  year  i (c^  ^,y^  j)  be  such  that 

the  90%  fill  rate  criterion  is  exceeded  (this  may  happen  if  population 
size  is  decreasing,  that  is,  units  are  being  retired  from  service),  then 
the  algorithm  immediately  goes  into  the  reduction  mode,  first  trying  to 
reduce,  to  the  greatest  extent  possible,  spares  (assuming  that  they  are 
more  expensive  than  servers  and  hence  should  have  a larger  salvage  value) 
and  then  attempting  to  reduce  servers  to  get  as  close  to  the  90%  boundary 
line  as  possible.  Again,  if  servers  are  more  expensive  than  spares,  the 
words  "server"  and  "spare"  should  be  interchanged. 

4.  Sample  Results 

Using  the  marine  gas  turbine  fleet  as  an  example.  Table  1 shows 
the  results  for  an  11-year  planning  horizon,  from  1975  to  1983,  using  an 
interest  rate  of  10%.  The  fleet  size  starts  out  at  a relatively  small 
number,  building  up  to  full  strength  by  year  1985.  The  component  exhibited 
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Figure  1 - Heuristic  algorithm. 

is  the  gas  generator  component.  Similar  calculations  were  performed  for 
the  power  turbine  component.  Generators  and  turbines  are  completely  in- 
dependent as  they  require  different  repair  facilities,  and,  of  course, 
spares  are  not  interchangeable. 

The  first  nine  columns  of  Table  1 (except  the  X column.  Column  4) 
are  input.  Column  1 gives  the  year,  Column  2 the  anticipated  population  size. 
Column  3 the  mean  failure  rate  schedule  for  the  particular  CIP  chosen  (in 
failures  per  day),  Column  5 the  average  service  (including  removal,  trans- 
portation, and  repair)  time  (in  days),  and  Columns  6-9  the  purchase  cost  of 
servers,  purchase  cost  of  spares,  unit  repair  cost,  and  investment  in  CIP 
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Column  4 and  the  last  five  columns  are  output.  It  was  convenient 
to  show  X^  (Column  4)  next  to  X^  so  that  one  can  readily  see  the  reli- 
ability growth.  Columns  10  and  11  give  the  algorithm's  "best"  c^,y^  • 

Column  12  shows  average  number  of  units  repaired  for  each  year,  while 
Column  13  provides  the  expenditures  (in  thousands  of  dollars)  up  to  and 
including  year  i.  Column  14  gives  the  present  worth  of  the  discounted 
cost  stream  up  to  and  including  year  i.  The  final  value  in  Column  14  is 
the  present  worth  of  all  expenditures  over  the  complete  planning  horizon, 
again  in  thousands  of  dollars. 

Another  run  is  presented  in  Table  2 which  differs  only  in  the  CIP 
chosen.  This  illustrates  how  one  can  study  the  consequence  of  various 
alternative  CIP's.  There  we  assume  that  in  order  to  save  CIP  investment 
cost  (C^),  we  will  stop  the  program  at  year  1978,  thus  achieving  only  a 

minimum  component  failure  rate  of  0.001  failures  per  day.  By  so  doing, 
we  save  the  two  large  expenditures  in  1979  and  1980,  and  Immediately  drop 
to  the  $350,000  per  year  needed  to  maintain  the  achieved  reliability  of 
0.001  failures  per  day.  However,  in  comparing  the  present  worth  of  the 
cost  streams  for  the  two  cases,  we  see  that  it  actually  costs  about  $5 
million  more  to  do  this  since  we  require  more  servers  and  spares.  Thus 
in  this  case  the  added  CIP  investment  is  well  worth  the  expenditure. 

5.  Program  Perturbation  Capability 

The  heuristic  algorithm  operates  on  a year-by-year  basis  and  hence 
has  the  limitation  that  it  does  not  "look  ahead."  While  we  believe  the 
heuristic  algorithm  does  a good  job  in  giving  an  optimal  or  near  optimal 
point  (c,y)  for  any  year,  when  considering  the  entire  planning  horizon 
it  may  not  prove  as  effective.  For  example,  consider  Table  1 once  again, 
particularly  year  1980.  We  see  for  that  year,  a jump  from  12  servers  in 
1979  to  15  in  1980  is  required,  but  in  the  following  year  only  13  servers 
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are  needed.  Recall  that  no  explicit  salvage  value  is  received  for  decreas- 
ing servers.  This,  In  Itself,  is  not  too  serious,  except  we  note  that  for 
1984  we  must  increase  servers  back  up  to  15,  having  to  purchase  two  addi- 
tional servers  for  that  year.  Had  we  kept  the  15  from  1981  we  would  not 
have  needed  to  make  further  server  purchases. 

But  more  importantly,  perhaps  we  could  have  avoided  going  to  15 
servers  in  1980  had  we  purchased  an  additional  spare,  especially  since  it 
is  needed  for  year  1981  anyway.  The  algorithm  allows  us  to  go  back  to  1980 
and  use  as  a starting  point  c=12  , y=13  (instead  of  the  1979  solution 
c=12  , y=ll  ).  This  may  give  a different  c,y  schedule  from  that  year  on, 
which  might  have  a cheaper  fixed  present  worth.  This  capability  to  go  back 
to  any  year  in  the  schedule  and  use  a starting  point  other  than  the  previous 
year's  c,y  we  refer  to  as  the  perturbation  capability.  Using  a different 
starting  point  (if  chosen  "intelligently")  might  get  us  to  a different  fea- 
sible solution  point  which  avoids  temporary  increases  in  either  c or  y . 
The  program  then  continues  its  calculation  from  the  year  for  which  the  per- 
turbation is  desired,  adding  the  new  cost  stream  with  the  correct  discount 
factor  from  that  year  forward,  giving  us  the  present  worth  for  the  new 
schedule  which  may  result.  Table  3 shows  this  effect  for  perturbing  the 
Table  1 solution  at  year  1980,  using  the  starting  point  c=12  , y=13  rather 
than  the  1979  solution  of  c=12  , y=ll  . We  see  for  the  new  resulting 
schedule  (which  differs  only  for  year  1980)  that  the  total  present  worth  is 
now  about  $39.05  million,  a savings  of  approximately  $0.1  million.  Further 
possible  perturbations  are  also  shown  with  the  resulting  costs.  No  other 
obvious  perturbations  are  indicated,  although  earlier  purchasing  to  avoid 
inflation  might  be  evaluated. 


At  first  glance,  it  appears  this  multi-year  programming  problem 
might  be  amenable  to  a dynamic  programming  solution  where  the  years  are 


stages  and  c and  y are  the  state  variables.  Aside  from  the  fact  that 
there  are  two  state  variables  (which  makes  the  computations  formidable) 
the  problem  is  not  decomposable.  The  fill  rate  constraint  must  hold  for 
each  year.  This  constraint  (for  year  i,  say)  Involves  the  q . which 
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decision  variables  for  year  i (c  ,y  ) , but  all  previous  c's  and  y's 

1 1 



since  A is  a function  of  A which  Is  a function  of  A.  „ , etc., 
i i—  I l~  Z 

[see  Equation  (2)].  This  prevents  starting  at  the  last  stage  and  pro- 

* 

ceeding  backward.  Only  if  the  reliability  were  constant  throughout  the 
multi-year  planning  horizon  could  dynamic  programming  be  attempted.  How- 
ever, we  believe  that  the  heuristic  algorithm,  coupled  with  the  perturba- 
tion capability,  does  provide  good,  although  not  necessarily  optimal, 
schedules. 


6.  Series  Queueing  Model 

The  queueing  model  portion  of  the  program  given  by  Equation  (1) 
assumes  that  removal,  transportation,  and  repair  are  a single  service 
action  so  that  if  all  servers  are  busy  and  a component  fails,  it  must 
wait  in  a queue  prior  to  removal  for  an  available  server.  Once  removal 
is  initiated,  no  further  queues  for  transportation  or  repair  are  encountered. 

In  many  situations,  the  removal,  transportation,  and  repair  phases 

may  be  separate,  each  having  their  own  "servers"  with  their  own  queues. 

To  model  this  situation  would  require  the  p . and  q . to  be  deter- 

n,i  n,i 

mined  by  a closed  network  or  cyclic  queueing  model  which  can  account  for 
the  finiteness  of  the  population,  but  requires  quite  involved  calculations 
(see,  for  example,  [3],  [10],  and  [11]).  The  question  then  naturally  arises 
as  to  how  crucial  the  finite  source  assumption  might  be,  since  rather  simple 
queueing  formulas  exist  for  infinite  source  series  queues  with  exponential 
arrival  and  service  patterns.  For  example,  suppose  for  component  removal, 
c^  removal  teams  are  available  for  the  system  (c^  could  equal  the  popu- 
lation size  N in  which  case  we  have  an  ample  server  model).  Assuming 
removal  times  are  exponential  and  failures  are  Poisson,  the  removal  portion 
can  be  modeled  by  an  M/M/c^  queue.  The  output  of  such  a queue  (assuming  an 

infinite  population)  is  Poisson.  Further,  suppose  there  were  c^  transport 

vehicles  available  to  the  system  for  shipping  components  to  the  repair  depot. 
The  transportation  phase  would  then  be  an  M/M/c?  queue,  assuming  transport 
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times  are  exponential.  Finally,  assuming  repair  channels  (again 

with  an  infinite  population),  the  repair  portion  becomes  an  M/M/c^  queue. 


Letting  p^  ^ be  the  probability  of  n components  at  the  jth 

phase  (j=l,2,3)  for  year  i,  it  can  be  shown  (see  Gross  and  Harris  [4], 
page  203)  that  the  joint  probability  of  1 in  the  removal  portion,  m 
in  the  transportation  portion,  and  n in  repair  for  year  i (call 


P^m.n;^  is  ®iven  ^ 


= n(1)  n (2)  n(3) 

p£,m,n;i  p H , i fm,i  pn,i 

Dropping  the  year  subscript  i for  the  time  being,  the  f il 1 rate  con- 
straint for  each  year  corresponding  to  Equation  (5)  becomes 


ry  r (i)  (2)  (3) 

2^  pa  pm  pn 


l m n 


> 0.90  , 


where  the  summation  is  taken  over  all  combinations  of  2,,m,n  such  that 

(2,+m+n  _<  y-1)  . Here  the  q£3  ^ equal  the  p£3  ^ , as  we  are  assuming  an 

infinite  source  model.  Since  p£3^  is  readily  computable  for  M/M/c 

models,  a series  queue  representation  is  possible  as  long  as  we  can  jus- 
tify approximating  the  finite  source  situation  by  an  infinite  source  model. 
Assuming  for  the  moment  we  can  (this  will  be  discussed  more  fully  below) , 
then  the  previous  methodology  can  be  used  if  we  substitute  for  the  pre- 
viously used  q . and  p . , the  sum  of  probabilities 
n,  l n, i 


St+m+k=n 


a)  (2)  (3) 

PJL,i  Pm,i  Pk,i  ’ 


with  the  arrival  rate  for  this  infinite  source  series  queueing  model  ad- 
justed to  be 

= NjTj.  . (1 


Once  again,  dropping  the  year  subscript  i for  convenience,  the 
standard  M/M/c  formulas  (see  Gross  and  Harris  [4J,  page  96,  for  example) 
give 
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with  X 
phase  j 


being  given  by  Equation  (7),  being  the  number  of  servers  for 

, and  Pj  being  the  mean  service  rate  for  phase  j (j=l  denotes 


removal,  j=2  denotes  transportation,  j=3  denotes  repair).  In  some  sit- 
uations, there  may  even  be  a fourth  phase  consisting  of  transportation  from 
repair  depot  to  spares  pool. 


The  mathematical  programming  problem  which  we  attacked  via  a heuristic 
algorithm  is  now  considerably  more  complex  if  c^  and  c^  are  also  decision 

variables,  since  we  now  seek,  for  each  year,  the  best  combination  (c^ , C2, 

c^,  y)  and  not  merely  (c,y)  . We  no  longer  have  a single  A , but  three 

A's  , say  A^,  A^,  A^  . For  equal  dollar  expenditures,  then,  we  can  buy  one 

spare  or  A^  removal  teams  or  A^  transport  vehicles  or  A^  repairmen. 

Further,  there  are  combinations  such  as  aAj  removal  teams  + 0A^  transport 

vehicles  + yA^  repairmen  (0  £ a,B,Y  £ 1)  , which  might  equal  the  purchase 

cost  of  a spare.  Conceptually,  the  extension  to  our  heuristic  algorithm 
would  be  that  given  a particular  c^,  c^,  c^  and  y , calculate  the  in- 
crease in  fill  rate  when  moving  to  all  points  of  equal  expenditure,  that 
is,  calculate  the  fill  rate  for  (c^  Cy  Cy  y+1)  , (Cj+A^ , Cy  c.3,  y)  , 

(Ci«  c2+A2’  c3’  y^  » (ci’  c2’  c3+A3’  y^  and  C2+^A2’  C3+Y^3’  y^ 

for  all  appropriate  ot,B,Y  • Cor  simplification,  the  equal  expenditure 
points  involving  combinations  of  c^,  c^,  c.^  (those  points  with  the  a. 
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3,  and  y terms)  might  be  ignored  so  that  at  each  step  we  add  either  one 
spare,  or  A^  removal  teams,  or  transports,  or  repairmen. 

If,  on  the  other  hand,  and  c are  not  decision  variables 

(and  in  our  applications  they  were  not),  then  the  algorithm  need  not  be 
modified  since  we  seek  only  "optimal  points"  (c^,y)  . Thus,  once  the 

p^'s  and  cln's  are  replaced  by  their  series  model  counterparts  as  given 

in  Equations  (6),  (8)  and  (9),  the  algorithm  is  exercised  as  before,  with 
c^  and  c-2  being  fixed  values.  In  the  next  section  we  present  some  re- 
sults for  such  a case  after  first  considering  the  question  of  accuracy  when 
using  an  infinite  source  approximation  to  a finite  population. 

7 . Accuracy  of  Infinite  Source  Approximations 

In  order  to  use  a series  queueing  model,  the  effect  of  assuming  an 
infinite  population  for  calculating  state  probabilities  when  in  actuality 
the  population  is  finite  must  be  investigated.  Let  us  consider  now  the 
non-series  finite  source  repairman  with  spares  model  given  by  Equation  (1) 
and  its  infinite  population  counterpart,  an  M/M/c.  model  witli  an  arrival 
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_ 00  y y+N 

\ = 23  \_Pn  = £ NXp  + 53  (N-n+y) Ap 

n=0  n n=0  n=y+l 

y+N  _ y+N 

= 23  N>lPn  " 21  (n-y)Xp 

n=0  n=y+l 

_ __  y+N  _ I"  y+N 

= nx  - x 53  (n-y)p„  = x n - 51  (n_y)pn 

n=y+l  n L n=y+l 

which  incidentally  is  how  R of  Equation  (3)  is  derived.  Thus  if 
y+N 

53  (n-y)p  is  small  in  relation  to  N , the  infinite  source  approxima- 
n=y+l  n 

tion  should  be  adequate  since  X = NX  , the  quantity  assumed  by  the  infi- 
nite source  model.  This  should  be  the  case  for  systems  with  not  too  much 

congestion  (p  small  for  n > y)  , y and  N large.  The  fact  that  a 

fill  rate  of  90%  must  be  guaranteed  should  require  c and  y large  enough 
such  that  there  is  a relatively  small  amount  of  congestion  in  the  system, 
thus  making  the  approximation  good  even  for  small  N . 

The  first  and  third  blocks  of  Table  4 show  comparisons  for  the  two 

cases  run  in  Tables  1 and  2,  respectively,  with  an  infinite  source  M/M/c 

used  in  calculating  p . . The  results  differ  only  slightly,  with  the 

n,  i 

infinite  source  model  requiring  an  additional  server  or  two  in  a few  of 

the  years.  The  infinite  source  model  will  occasionally  require  slightly 

higher  c (or  y)  than  the  finite  source  model  in  order  to  meet  the  avail- 
ability constraint.  This  is  caused  by  the  lack  of  state  dependence  in  the 
arrival  rate  for  the  infinite  source  model  which  has  the  effect  of  making 
the  infinite  source  arrival  rate  slightly  higher  than  the  arrival  rate  for 
the  corresponding  finite  source  model.  Hence  the  infinite  source  model  is 
a conservative  approximation  with  respect  to  fill  rate. 

Since  the  infinite  source  approximation  appears  to  give  good  re- 
sults, we  have  run  an  infinite  source  three-stage  (removal,  transport, 
and  repair)  series  model,  using  the  same  data  as  in  Tables  1 and  2,  but 
assuming  removal  is  M/M/°°,  transport  is  M/M/00,  and  repair  is  M/M/c  (c  to 
be  determined  as  before). 
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Actually,  for  ample  server  queues.  It  is  not  necessary  to  assume 
exponential  service  since  the  output  of  M/G/°°  is  also  Poisson  (see,  for 
example,  Mirasol  [7]).  Thus,  these  results  also  hold  for  M/G/°°  removal 
and  transport  phases.  The  mean  single  stage  repair  time  is  allocated  to 
the  multi-stage  model  so  that  on  the  average  5%  is  consumed  in  removal, 

20%  in  transportation,  and  75%  in  repair.  These  results  comparing  the 
two  models  are  also  shown  in  Table  4.  Here,  in  many  of  the  years,  fewer 
servers  are  required  to  achieve  the  90%  fill  rate  when  treating  removal, 
transportation,  and  repair  separately.  This  is  as  we  would  expect,  since 
if  these  three  phases  are  separate,  removal  and  transportation  can  occur 
even  if  there  is  a queue  at  the  repair  facility.  By  treating  the  separate 
phases  together  as  one  service  operation,  as  in  the  single-stage  model,  we 
are  in  essence  requiring  a free  repair  facility  prior  to  removal , in  which 
case  the  repairman  would  be  idle  until  the  component  arrived.  Thus  in  the 
situation  where  removal  and  transportation  have  ample  servers  and  queueing 
takes  place  only  at  repair,  conservative  results  are  obtained  if  we  treat 
the  three  separate  phases  as  one.  If  indeed  the  service  portion  is  made  up 
of  separate  phases,  a more  realistic  representation  of  the  system  is  obtained 
by  using  the  series  model  in  spite  of  the  slight  inaccuracy  caused  by  in- 
voking the  infinite  source  assumption  when  calculating  probabilities. 

With  fill  rate  constraints  of  90%  or  better  which  guarantee  relatively 
little  congestion,  it  appears  that  the  inaccuracies  caused  by  using  the 
infinite  source  assumption  in  the  series  model  are  negligible  as  compared 
to  those  that  result  when  using  a single-stage  finite  source  model  if  the 
service  actually  is  multi-stage. 
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